REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 

of  information,  including  suggestions  for  reducing  this  burden  to  Washington  Headquarters  Service,  Directorate  for  Information  Operations  and  Reports, 

1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget, 

Paperwork  Reduction  Project  (0704-0188)  Washington,  DC  20503. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


5b.  GRANT  NUMBER 

FA9550-08-1-0173 

5c.  PROGRAM  ELEMENT  NUMBER 


5e.  TASK  NUMBER 

N/A 

5f.  WORK  UNIT  NUMBER 

N/A 


12.  DISTRIBUTION  AVAILABILITY  STATEMENT 

Distribution  A 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  research  performed  during  this  AFOSR  grant  has  extended  the  basic  discrete  DMOC  (Discrete  Mechanics  and 
Optimal  Control)  framework  to  derive  algorithms  for  general  nonholonomic  multi-body  systems,  with  applications  to  a 
variety  of  problems  including  space  mission  design.  A  computational  framework  was  developed  which  can  automatically 
construct  integration  and  optimization  schemes  by  providing  a  high-level  description  of  the  mechanical  system  in  terms 
of  its  physical  layout,  inertial  properties,  constraints,  actuation,  and  external  influences.  This  last  year  of  the  grant  was 
focused  on  optimization  problems  with  non-holonomic  constraints  as  well  as  continued  applications  to  space  mission 
design.  Besides  describing  the  latest  results  and  publications  obtained  over  the  past  12  months,  we  also  describe  the 
next  steps  that  we  expect  to  explore  now  that  the  work  of  this  grant  is  over. 

15.  SUBJECT  TERMS 


17.  LIMITATION  OF  18.  NUMBER 
ABSTRACT  OF  PAGES 

SAR  3 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Mathieu  Desbrun 

19b.  TELEPONE  NUMBER  (Include  area  code) 

626  395  6230 


16.  SECURITY  CLASSIFICATION  OF: 
a.  REPORT  I  b.  ABSTRACT  I  c.  THIS  PAGE 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

California  Institute  of  Technology 
1200  E.  California  Blvd 
Pasadena  CA  91125 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Office  of  Scientific  Research 
875  N.  Randolph  St 
Rm  3112 

Arlington,  VA  22203 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

N/A 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

AFOSR 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

AFRL-OSR- VA-TR-20 1 2-0867 


5d.  PROJECT  NUMBER 

N/A 


3.  DATES  COVERED  (From  -  To) 

04/01/2008-  11/30/2010 

5a.  CONTRACT  NUMBER 


2.  REPORT  TYPE 

Final  Technical  Report 

4.  TITLE  AND  SUBTITLE 

Geometric  Computational  Mechanics  and  Optimal  Control 


6.  AUTHOR(S) 

Mathieu  Desbrun,  Marin  Kobilarov 


1.  REPORT  DATE  (DD-MM-YYYY) 

12/02/2011 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI-Std  Z39-18 


INSTRUCTIONS  FOR  COMPLETING  SF  298 


1.  REPORT  DATE.  Full  publication  date, 
including  day,  month,  if  available.  Must  cite  at 
lest  the  year  and  be  Year  2000  compliant,  e.g., 
30-06-1998;  xx-08-1 998;  xx-xx-1 998. 

2.  REPORT  TYPE.  State  the  type  of  report,  such 
as  final,  technical,  interim,  memorandum, 
master's  thesis,  progress,  quarterly,  research, 
special,  group  study,  etc. 

3.  DATES  COVERED.  Indicate  the  time  during 
which  the  work  was  performed  and  the  report 
was  written,  e.g.,  Jun  1997  -  Jun  1998;  1-10  Jun 
1996;  May  -  Nov  1998;  Nov  1998. 

4.  TITLE.  Enter  title  and  subtitle  with  volume 
number  and  part  number,  if  applicable.  On 
classified  documents,  enter  the  title  classification 
in  parentheses. 

5a.  CONTRACT  NUMBER.  Enter  all  contract 
numbers  as  they  appear  in  the  report,  e.g. 

F3361 5-86-C-51 69. 

5b.  GRANT  NUMBER.  Enter  all  grant  numbers 
as  they  appear  in  the  report,  e.g. 

1 F665702D1257. 

5c.  PROGRAM  ELEMENT  NUMBER.  Enter  all 
program  element  numbers  as  they  appear  in  the 
report,  e.g.  AFOSR-82-1234. 

5d.  PROJECT  NUMBER.  Enter  al  project 
numbers  as  they  appear  in  the  report,  e.g. 

1 F665702D1257;  ILIR. 

5e.  TASK  NUMBER.  Enter  all  task  numbers  as 
they  appear  in  the  report,  e.g.  05;  RF0330201; 
T4112. 

5f.  WORK  UNIT  NUMBER.  Enter  all  work  unit 
numbers  as  they  appear  in  the  report,  e.g.  001; 
AFAPL30480105. 

6.  AUTHOR(S).  Enter  name(s)  of  person(s) 
responsible  for  writing  the  report,  performing  the 
research,  or  credited  with  the  content  of  the 
report.  The  form  of  entry  is  the  last  name,  first 
name,  middle  initial,  and  additional  qualifiers 
separated  by  commas,  e.g.  Smith,  Richard,  Jr. 

7.  PERFORMING  ORGANIZATION  NAME(S) 
AND  ADDRESS(ES).  Self-explanatory. 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER.  Enter  all  unique  alphanumeric  report 
numbers  assigned  by  the  performing 
organization,  e.g.  BRL-1234;  AFWL-TR-85- 
401 7-Vol-2 1  -PT -2 . 

9.  SPONSORING/MONITORS  AGENCY 
NAME(S)  AND  ADDRESS(ES).  Enter  the  name 
and  address  of  the  organization(s)  financially 
responsible  for  and  monitoring  the  work. 

10.  SPONSOR/MONITOR'S  ACRONYM(S). 

Enter,  if  available,  e.g.  BRL,  ARDEC,  NADC. 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S).  Enter  report  number  as  assigned 
by  the  sponsoring/  monitoring  agency,  if 
available,  e.g.  BRL-TR-829;  -215. 

12.  DISTRIBUTION/AVAILABILITY 
STATEMENT.  Use  agency-mandated 
availability  statements  to  indicate  the  public 
availability  or  distribution  limitations  of  the  report. 
If  additional  limitations/restrictions  or  special 
markings  are  indicated,  follow  agency 
authorization  procedures,  e.g.  RD/FRD, 

PROPIN,  ITAR,  etc.  Include  copyright 
information. 

13.  SUPPLEMENTARY  NOTES.  Enter 
information  not  included  elsewhere  such  as: 
prepared  in  cooperation  with;  translation  of; 
report  supersedes;  old  edition  number,  etc. 

14.  ABSTRACT.  A  brief  (approximately  200 
words)  factual  summary  of  the  most  significant 
information. 

15.  SUBJECT  TERMS.  Key  words  or  phrases 
identifying  major  concepts  in  the  report. 

16.  SECURITY  CLASSIFICATION.  Enter 
security  classification  in  accordance  with 
security  classification  regulations,  e.g.  U,  C,  S, 
etc.  If  this  form  contains  classified  information, 
stamp  classification  level  on  the  top  and  bottom 
of  this  page. 

17.  LIMITATION  OF  ABSTRACT.  This  block 
must  be  completed  to  assign  a  distribution 
limitation  to  the  abstract.  Enter  UU  (Unclassified 
Unlimited)  or  SAR  (Same  as  Report).  An  entry  in 
this  block  is  necessary  if  the  abstract  is  to  be 
limited. 


STANDARD  FORM  298  Back  (Rev.  8/98) 


Annual  Performance  Report,  2011 

Grant  No:  FA9550-08-1-0173 
Geometric  Computational  Mechanics  and  Optimal  Control 
(Pis:  Marsden/Desbrun) 

The  last  year  of  this  grant  continued  the  development  and  numerical 
implementation  of  DMOC  (Discrete  Mechanics  and  Optimal  Control) 
with  applications  to  a  variety  of  optimization  problems.  Professor  J.E. 
Marsden,  PI  on  this  grant,  unfortunately  passed  away  this  past 
September;  Prof.  Desbrun,  who  has  collaborated  with  Marsden  and 
his  postdocs  for  the  past  six  years  on  DMOC-related  projects,  has 
taken  over  the  task  of  completing  this  research  effort. 

This  year  was  focused  on  optimization  problems  with  non-holonomic 
constraints  as  well  as  continued  applications  to  space  mission  design. 
The  method  builds  on  what  is  by  now  well  developed  and  successful 
techniques  for  the  theoretical  and  numerical  implementation  of 
discrete  mechanics  done  over  the  last  decade.  DMOC  combines 
these  discrete  mechanics  techniques  with  direct  methods  for  optimal 
control,  such  as  SQP  (sequential  quadratic  programming)  and  related 
methods.  Further  methods  that  depend  on  global  optimization 
problems  are  in  development  and  preliminary  versions  of  these 
results,  many  of  which  were  obtained  with  the  current  postdoctoral 
fellow  Marin  Kobilarov,  were  reported  in  the  recent  annual  meeting. 

This  grant’s  achievements  were  done  not  only  with  graduate  students,  but 
also  with  the  recent  past  AFOSR  (partially)  supported  postdoc,  Sigrid 
Leyendecker,  a  recent  winner  of  an  Emmy  Noether  5  year  research  award 
to  start  her  group  in  Germany  and  Sina  Ober-Blobaum,  now  a  junior 
professor  at  Paderborn.  Our  collaborations  on  DMOC  are  continuing,  for 
example  with  Ober-Blobaum  on  multiscale  and  receding  horizon  versions 
of  DMOC  with  applications  to  circuits,  and  with  Leyendecker  on  3D 
multibody  simulations  and  design  of  dynamics. 

The  concrete  progress  that  was  made  since  the  last  reporting  period 
include: 


Basic  Optimal  Control  Theory  and  Applications 

This  line  of  work  is  based  on  variational  methods  for  time  evolution  for 
constrained  systems,  which  showed  a  marked  improvement  over  existing 
methods.  The  techniques  devised  in  previous  years  of  this  grant  have  also 
been  extended,  both  time  evolution  and  optimal  control  problems,  to  non- 
holonomic  constraints  and  Lie  groups: 

M.  Kobilarov,  J.  Marsden,  Discrete  Geometric  Optimal  Control  on  Lie  Groups , 
IEEE  Transactions  on  Robotics  and  Automation,  to  appear,  2011. 

M.  Kobilarov,  D.  Martin  de  Diego,  S.  Ferraro,  Simulating  Nonholonomic 
Dynamics.  Boletin  de  la  Sociedad  Espanola  de  Matimatica  Aplicada  (SeMA),  50, 

2010,  pp  61-81. 

K.  FlaBkamp,  S.  Ober-Blobaum,  M.  Kobilarov,  Solving  optimal  control  problems 
by  using  inherent  dynamical  properties,  Proceedings  in  Applied  Mathematics  and 
Mechanics,  2010,  vol.  10-1,  pp.  577-578 

S.  Dubljevic,  M.  Kobilarov,  and  James  Ng,  Discrete  Mechanics  Optimal  Control 
(DMOC)  and  Model  Predictive  Control  (MPC)  Synthesis  for  Reaction-Diffusion 
Process  System  with  Moving  Actuator,  American  Control  Conference  (ACC), 

2010,  pp.  5694-5701. 

Applications  to  Space  Mission  Design. 

The  optimization  of  spacecraft  trajectories  has  been  a  natural  application 
area  for  DMOC.  We  published  a  paper  this  year  building  upon  last  year’s 
spacecraft  trajectory  optimization  paper  with  DMOC  to  include  refinement 
strategies  for  faster  convergence: 

A.  Moore,  Sina  Ober-Blobaum,  and  J.E.  Marsden,  Mesh  Refinement  Strategies  for 
Spacecraft  Trajectory  Optimization  using  Discrete  Mechanics  and  Optimal 
Control,  AAS  1 1-127,  2011 


Conclusions  &  Next  Steps 

The  research  performed  during  this  AFOSR  grant  extended  the  basic 
discrete  DMOC  framework  to  derive  algorithms  for  general  nonholonomic 


multi-body  systems.  A  computational  framework  was  developed  which  can 
automatically  construct  integration  and  optimization  schemes  by  providing 
a  high-level  description  of  the  mechanical  system  in  terms  of  its  physical 
layout,  inertial  properties,  constraints,  actuation,  and  external  influences. 
The  main  focus  was  then  to  extract  and  factor  out  as  much  structure  as 
possible  to  obtain  algorithms  with  lowest  dimension  that  are  amenable  to 
controllability  and  convergence  analysis,  i.e.,  determining  whether  the 
system  can  reach  a  desired  state  and  achieve  it  in  an  optimal  way. 

In  this  context  a  novel  optimal  control  method  was  developed  to  exploit  the 
structure  of  optimal  trajectories.  As  a  result,  instead  of  relying  on  a  black¬ 
box  optimization  package,  optimal  control  trajectories  are  derived  through  a 
discrete  optimality  principle  similar  to  the  construction  of  higher  order 
geodesics  in  the  continuous  setting. 

Consequently,  globally  optimal  methods  for  computing  optimal  trajectories 
for  vehicles  with  complex  dynamics  were  developed.  The  basic  idea  is  to 
use  motions  optimized  using  discrete  optimal  control  in  order  to  generate 
near-optimal  trajectories  which  probabilistically  explore  the  state  space 
searching  for  the  optimal  path.  The  set  of  all  such  trajectories  are 
approximated  using  a  graph  with  nodes  corresponding  to  randomly 
sampled  states.  The  methodology  is  designed  to  converge  to  an  optimal 
solution  as  the  number  of  nodes  is  increased.  The  approach  was 
successfully  applied  to  a  variety  of  systems  in  simulation  as  well  as  to  a 
real  vehicle  --  an  unmanned  helicopter  flying  in  an  urban  terrain.  This 
promising  research  direction  will  now  be  further  studied  in  the  context  of 
future  grants. 


